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^ ■. Abstract 

f~^ ' We show that FFE (Force-Free Electrodynamics) describes pulsar magnetospheres - even when 

(N ■ 

j_^ ■ radiation and particle wind power are comparable to the spin-down power. At the same time, we 

^ , show that FFE is insufficient for calculating pulsar magnetospheres. This is because FFE admits 

C*^ I a large family of stationary solutions with different currents flowing in the closed line region. To 
choose the actual solution one needs a model of current generation which goes beyond pure FFE. 

^ , We calculate several FFE magnetospheres for the aligned rotator. The Poynting power of these 

^^ \ solutions fills the range c~^/i^O^ < L < 0.67c^^ fi^^^{c/^Rs)'^ , for angular velocity O and 

f~^ ' magnetic field which is a pure dipole /i on the surface of the star of radius Rg. 

^^ . Anomalous braking indices of young pulsars might be explained by currents flowing in the closed- 
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line region. 



I. INTRODUCTION 



It might appear that pulsar magnetosphere and spin-down power have fina 
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ated. Both aligned rotator (pulsar with small spin-dipo 
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9|, |10[ and inclined rotator (arbitrary spin-dipole angle) 
magnetospheres were calculated in the framework of force-free electrodynamics (FFE, §2). 
Magnetohydrodynamics (MHD) codes were also used, but we will argue that MHD applies 
to pulsars only inasmuch as it reduces to FFE. 

While we do think that FFE is applicable to pulsars (§2), we also show that FFE de- 
scription of pulsars is ambiguous. In particular, we show (§3) that FFE can only give lower 
and upper bounds for the pulsar Poynting power: 

~ -^ ~ „3 o D / ' y^J 



where Rs is the radius of the star. 

Observations might actually favor replacing the standard power formula, L ~ c~^/i^fi^, 
with the bounds Q. Anomalous braking indices of young pulsars are best explained by the 
spin-down power L oc i7^~" [1^. This spin-down power obtains if we assume that Poynting 
and spin-down power are comparable, and also assume that L is proportional to some power 
of Q which falls into the range allowed by FFE bounds. Then we get the spin-down power 
L ~ c^^fj,'^Q^{c/QRs)°', with < a < 2, in agreement with observations. 

II. FORCE-FREE ELECTRODYNAMICS OF PULSARS 

Here we show that FFE describes large-scale electromagnetic field of pulsars, even when 
particle or radiation power are comparable to the large-scale Poynting power. 
FFE postulates that electromagnetic field satisfies Maxwell equations, 

dtB = -VxE, dtE = VxB-j, (2) 

with the current 
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(B-VxB-E-Vx E)B + (V ■ E)E x B 

Initial condition E ■ B = is assumed, and then this condition is sustained by FFE. 

The strange-looking Ohm's law (jH)) is designed [l^ in such a way that force-free condition 



pE + j X B = (4) 



(with p = V • E) is always satisfied. We will explain physical meaning of the FFE Ohm's 
law below. 

Standard derivation of FFE postulates that particle inertia is negligible compared to 
characteristic values of electromagnetic forces. Then the electromagnetic force density should 
vanish, as given by the force-free condition (j^. One then finds an expression for the current 
j which sustains this condition as the fields evolve according to Maxwell equations (J2I). 

The force-free condition is Lorentz- invariant: F^^j^ = 0, where F'^'^ is the electromag- 
netic field tensor and j^ is the 4-current. Since divergence of the electromagnetic energy- 
momentum tensor is d^T^^ = —F^'^j^, FFE should apply if energy-momentum tensor is 
dominated by the electromagnetic part. 

But this condition - negligible energy-momentum tensor of particles - might be too 
.estnct.ve. Fo. example. X-.ay lun,ino=ity of C.ab-.ike pulsa. ca. be up to 10% of =pm- 
down power 14]. Since X-ray energy is coming from particles, the particle power might be 
actually comparable or even dominate the spin-down power. If so, the derivation of FFE 
presented above fails. However, we think that even in such cases FFE actually remains valid, 
because FFE equations can be derived without postulating small particle inertia. 

We start with a definition. Consider arbitrary electromagnetic field. At each event (space- 
time event), one can find a good reference frame, where the electric field E is parallel to the 
magnetic field B. In fact, there is a whole one-parameter family of good frames, obtained 
by boosting a good frame along common direction of E and B. Which good frame is chosen 
is irrelevant for our purposes. 

FFE follows from two assumption. In a good frame: (i) electric field vanishes, (ii) current 
flows along B. These conditions obviously give force freedom (0]) in a good frame, and 
therefore in any frame, and thus these conditions are equivalent to FFE. The physical 
meaning of FFE Ohm's law is most clear in a good frame: the current flows along B, 
and the magnitude of the current is chosen in such way as to sustain zero electric field. We 
will now argue that these conditions should indeed be satisfied, at least in Crab-hke and 
Vela-like pulsars. 

The first condition means that electric field component along magnetic field at radius r 
should be much smaller than the characteristic value ~ {Qr/c)B. Consider young pulsar 
with nominal surface field 5^ ~ 3 x 10^^ G and angular velocity Q ~ 100 s"^. Then 
electron-positron pair creation avalanche (|l^: electrons and positrons produce gamma-rays 



by curvature radiation, gamma rays produce pairs on magnetic field) discharge the parallel 
component of electric field at distances ^ 600km. At larger distances magnetic field becomes 
too weak for pair production. The light cylinder radius is larger, c/Vt = 3000km, and one 
might worry about shortening the field near the light cylinder. It is thought, however, that 
avalanches are producing dense outflows of electron-positron plasma [15[ . With the electron- 
positron plasma density well above Goldreich- Julian density, the parallel component of the 
electric field should shorten near the light cylinder too. 

The second condition - in a good frame current flows along magnetic field - should be 
easily satisfied. Fast synchrotron/cyclotron cooling puts charged particles into the lowest 
Landau level (cooling is fast not only for electrons and positrons, but also for ions, if there 
are any ions in the magnetosphere) . The magnetic dipole moment corresponding to the 
lowest Landau level is small and the resulting large scale perpendicular current is negligible. 

Thus: negligible particle inertia is not a must for FFE applicability. It is sufficient to 
have (i) easily available charges (to kill E along B), (ii) strong enough magnetic field (to 
ensure fast perpendicular cooling). It appears that these conditions can be satisfied, at least 
in Vela-like and Crab-like pulsars. 

III. STATIONARY FFE MAGNETOSPHERE AND POYNTING POWER 

Stationary FFE ma gne tosphere (electromagnetic field corotates with the star) satisfies 
the following equation [ll| : 

B X V X (B + Vx (Vx B)) = 0, (5) 

with V = fi X r and V • B = 0. The electric field is E = -V x B. 

Solutions of the pulsar magnetosphere equation © follow from the isotopological varia- 
tion variation principle: 6S = 0, where S is proportional to electromagnetic action 

S= fd\{B^-(VxBy), (6) 

and the variation of action should be calculated for isotopological perturbation of magnetic 
field 

6B = V x{6^xB). (7) 

In other words, action should be extremized by continuously displacing magnetic field lines. 
The footpoints of the lines remain fixed. The number density of footpoints is given by the 



normal component of magnetic field on the surface of the star Bn{6, 0), where 6 and are 
spherical angles. 

There are magnetic field lines of two types in the stationary magnetosphere: star-infinity 
and star-star. The star-star lines define a mapping of the closed-line part of the stellar 
surface onto itself. 

Stationary FFE magnetosphere can be calculated for a given (i) O, (ii) Bn{0,(f)), (iii) 
mapping of the closed-line part of the surface. The surface mapping should be consistent 
with Bn{6,(f)). But for a given Bn{6,(f)), there are infinitely many consistent mappings. 
Therefore there are infinitely many FFE magnetospheres with given ft and Bn{0,(f)). FFE 
has nothing to say about the actual mapping of the closed-line part of the stellar surface. 

Physically, different mappings correspond to different currents flowing in the closed line 
region. For axisymmetric pulsar, one usually postulates zero poloidal current in the closed- 
line region, but we see no reason why the microscopic mechanism that realizes FFE should 
obey this condition for a generic asymmetrical Bn{0). The small residual parallel electric 
field is adjusted by the requirement that charge density be equal to Goldreich- Julian. There 
is no reason to expect that this adjustment process terminates at zero current. Note also that 
FFE generically predicts singular current layers along magnetic separatrix (surface between 
open and close field lines). Again there is no reason to expect that microscopic physics 





FIG. 1: Standard axisymmetric pulsar. No current in the closed region. Separatrix at tpQ = 1.32. 
Power L = 1.08. (1) ^/i-isolines, with -(/'-intervals of O.l^o- (2) Poloidal current I{ip)- 



producing very high current densities near magnetic separatrix will not result in some finite 
current densities within the separatrix. 

We studied the effect of the inner currents numerically. We assumed axisymmetry. For 
stationary axisymmetric force-free magnetic field, the magnetosphere equation (jH)) reduces 
to the basic equation of [2] and gives magnetic field of the form 



B 
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(8) 



Here r, z are cylindrical coordinates, il){r,z)/r is the toroidal component of the vector 
potential, /(■?/') is the poloidal current, and ■?/'(r, z) satisfies the following equation 



1 2 

(1 - r^){:iprr + -'Ipr + ^^zz) - 'A + F{%Ij) = 0, 



(9) 



F{i;) = AIj^. The units are /i = c = fi = 1. 

To simplify the calculation, we assumed that negligibly small north-south asymmetries 
launch large closed currents. The numerical procedure is just the one discovered by j^. The 
only difference, is that now a non-zero toroidal field (corresponding to poloidal current) exists 
in the closed line region. Poloidal current in the open-line region is numerically adjusted 
to allow smooth crossing of the light cylinder. Poloidal current in the closed-line region is 
chosen in an arbitrary manner. 

The results are shown in figures (l)-(4). The radius of the star is Rg = 0.1. We have 
checked that we obtain legitimate FFE solutions with B > E. Figure (1) just repeats 




FIG. 2: With current in the closed region. Separatrix at ipo = 1.67. Power L = 1.75. 




FIG. 3: With current in the closed region. Separatrix at V'o = 2.36. Power L = 3.65. 

the standard solution. Figures (2) and (3) show solutions with increasing toroidal field in 
the closed region. As toroidal field increases, more lines open up and the pulsar power 
increases accordingly. With very large toroidal field in the closed region, one should get the 
split-monopole configuration shown in figure (4). 

We must acknowledge that our numerical procedure for calculating stationary magneto- 
spheres breaks down for large currents in the closed line region, and we cannot demonstrate 
the field configurations with V^o ~ 3. However we think that such solutions do exist. The 
action © is positive definite inside the light cylinder and therefore it must possess a mini- 
mum for arbitrarily large twist of the initial mapping. As the lines open up, the functional 
(jHl) gets decreased. 

A better way to argue that arbitrarily large opening of field lines can be achieved for large 
twists, is to consider a non-rotating FFE magnetosphere. Then equation Q simplifies. In 
spherical coordinates r = _Rsin6', z = Rcos9, 



^lJJ,J, + sine (^iJe) + F(^) = 0, 



(10) 



with selfsimilar solutions of the form 



V' oc /?-", / oc V 
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We show two such solutions in Fig. (5). Note, that the twist, defined as longitude difference 
between footpoints of the same field line, is not too big. Assumingly, if such a star is spun 




FIG. 4: All lines open. Separatrix at -00 = 10- Power L = 67.1 

up, it will generate Poynting flux proportional to magnetic energy density near the light 
cylinder, L oc i?^^^^""). 

We must conclude that FFE-allowed magnetospheres span a whole range between the 
standard magnetosphere with Poynting luminosity L ~ c~^/i^f2^, and split-dipole magneto- 
sphere with luminosity L ~ c~^/i^fi^(c/r2-Rs)^- Pure FFE cannot say which magnetosphere is 
actually realized. Further progress requires a model of current generation. Poloidal currents 
in the closed region might explain anomalous braking indices. 
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FIG. 5: ^-isolines of non-rotating star with poloidal current: (1) current I (x ip^, giving tp cc R ^' ^, 
maximal twist 118° (2) current / ex tp^, giving ip oc R~^'^, maximal twist 147° 



